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Abstract 

It has been considered that the effective bending rigidity of fluid membranes should be reduced 
by thermal undulations. However, recent thorough investigation by Pinnow and Helfrich revealed 
significance of measure factors for the partition sum. Accepting the local curvature as a statistical 
measure, they found that fluid membranes are stiffened macroscopically. In order to examine 
this remarkable idea, we performed extensive ab initio simulations for a fluid membrane. We set 
up transfer matrix which is diagonalized by means of the density-matrix renormalization group. 
Our method has an advantage in that it allows us to survey various statistical measures. As a 
consequence, we found that the effective bending rigidity flows toward strong coupling under the 
choice of local curvature as a statistical measure. On the contrary, for other measures such as 
normal displacement and tilt angle, we found clear tendency toward softening. 

PACS numbers: 87.16.Dg Membranes, bilayers, and vesicles, 87. 16. Ac Theory and modeling; computer 
simulation, 05.10.-a Computational methods in statistical physics and nonlinear dynamics, 05.10.Cc Renor- 
malization group methods 
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I. INTRODUCTION 



Amphiphilic molecules in water segregate spontaneously into flexible extended surfaces 
called fluid (lipid) membranes |l|, Q . The fluid membranes are free from both surface tension 
and shear modulus, and the elasticity is governed only by bending rigidity |3|, Jj|. The 
Hamiltonian is given by the following form, 

H = JdA (~ J 2 + kK\ . (1) 

The mean curvature J is given by the summation of two principal curvatures J = c\ + C2, 
whereas the Gaussian curvature K is given by their product K = c\Ci. The corresponding 
two moduli k and R are called bending rigidity and Gaussian-curvature modulus, respec- 
tively. The integration / dA extends over the whole membrane surface. Hereafter, we drop 
the Gaussian- modulus term (k = 0), because this term is topologically invariant |], §, and 
we restrict ourselves to a fixed topology (planar surface). 

Contrary to its seemingly simple form, the Hamiltonian (|l|) brings about perplexing 
problems: As a matter of fact, under an actual representation with differential geometry 
(for instance, see Eq. (|3]) explained afterwards), it turns out that membrane's undulations 
are subjected to complicated mutual scatterings. Therefore, it is expected that the effective 
bending rigidity for macroscopic length scales differs from the bare rigidity owing to the 
thermally activated undulations. Aiming to clarify this issue, a number of renormalization- 
group analyses have been reported so far [||, || |7j. The results are summarized in the 
following renormalizat ion-group equation, 

k' = k — a ——— log M (2) 

with renormalized bending rigidity k', temperature T and number of decimated molecules 
M. All literatures agree that the numerical prefactor in the above equation is a = 3. (More 
detailed account of historical overview would be found in Ref. ||.) Because of a > 0, 
effective bending rigidity is reduced by thermally activated undulations. This conclusion 
might be convincing, because the membrane shape itself should be deformed by thermal 
undulations. As a matter of fact, it has been known that the orientational correlation is lost 
for long distances M. It is quite natural to anticipate that membranes become floppy for 
those length scales exceeding this correlation length. 
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Recently, however, Pinnow and Helfrich || |10| obtained a remarkable conclusion a = 
— 1(< 0). The key ingredient of their new argument is that they considered the role of mea- 
sure factors for the partition sum. They insist that the local mean curvature J should be the 
right statistical measure rather than other measures such as normal displacement h(x, y) and 
local tilt angle 8(x,y). (The variable h(x,y) denotes the normal displacement of membrane 
from a base (reference) plane. We will explain it in the next section.) After elaborated cal- 
culation of the variable replacement h(x, y) — ► J(x, y) and succeeding renormalization-group 
analysis, the authors reach the conclusion of a = — 1. 

In order to examine their remarkable scenario, we performed first-principle simulation 
for a fluid membrane. Our simulation method has an advantage in that we cover various 
integration measures in a unified way. In addition, our simulation does not rely on any 
perturbative treatments. Hence, it would be meaningful to complement the analytical per- 
turbative treatment. We employed a scheme of transfer matrix, which is diagonalized by 



means of the density- matrix renormalization group |TT|, [L2|, [13], |14[ . It is to be noted that 
recently, such elastic (bosonic) systems come under through investigation by means of diag- 
onalization after the advent of density- matrix renormalization group [15, [1(| [17], [TBI [TjJ 20 



It has to be mentioned that the Monte-Carlo method has been utilized successfully in the 
studies of membranes and vesicles [2]. For the Monte-Carlo method, however, tethered 
(polymerized) membrane rather than fluid membrane is more suited EI| , because membrane 
is implemented in computer as an assembly of molecules and junctions bearing close resem- 
blance to tethered membrane. However, Gompper and Kroll succeeded in simulating fluid 
membranes by Monte-Carlo method, allowing reconstructions of junctions during simulation 
(dynamical triangulation) |22[] . They succeeded in observing the variation of topological in- 
dex with respect to the change of membrane concentration and temperature. In fairness, it 
has to be mentioned that their Monte-Carlo data indicate softening for lipid vesicles. 

The rest of this paper is organized as follows. In the next section, we explain our transfer- 
matrix formalism for a fluid membrane. We also explicate the density-matrix renormalization 
group with which we diagonalized the transfer matrix. In Section |TJ, we calculate the 
effective bending rigidity. For that purpose, we introduce a scheme of coarse-graining. Our 
data clearly support membrane stiffening under the choice of mean curvature J(x, y) as 
a statistical measure. For other measures, on the contrary, we observed tendency toward 
softening. In the last section, we present summary and discussions. 
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II. TRANSFER-MATRIX FORMALISM AND ITS DIAGONALIZATION 

THROUGH DENSITY-MATRIX RENORMALIZATION GROUP 



In this section, we present our numerical simulation technique. First, we explain our 
transfer-matrix formalism. Secondly, we introduce the density-matrix renormalization group 
with which we diagonalize the transfer matrix. Demonstration of the algorithm is also 
presented. 

A. Transfer-matrix formalism 

In order to describe the shape of membranes, it is convenient to use the Monge gauge ||23| . 
In this frame, membrane deformation is described by a normal displacement (deformation) 
h(x, y) from a base (reference) plane parameterized by Cartesian coordinates (x, y). In terms 
of this representation frame, the mean curvature J(x, y) is given by, 

_ (gft + 9 2 y h)(l + {8 X hf + (dyhf) - 2d X hdyhd X dyh ~ dfjh{d x hf ~ ^hjdyh) 2 

{x,y) (i + (d x hy + (d y hy)v* 

(3) 

explicitly. Similarly, the infinitesimal area dA is given by, 

dA = (1 + (d x h) 2 + {d y h) 2 ) 1/2 dxdy. (4) 

Putting all together into the Hamiltonian ([I]), we obtain an explicit expression in term of 
the displacement field h(x,y). 

Now, we are led to a two-dimensional scalar-field theory h(x,y). However, the theory is 
afflicted by very complicated interactions. The aim of this paper is to investigate the theory 
beyond perturbative level by means of an ah initio method: For that purpose, we put the 
theory on a square lattice with lattice constant a; see Fig. |l|. Accordingly, the field variables 
{hij} are now indexed by integer indices i and j. Throughout this paper, we set the lattice 
constant as the unit of length a = 1. In other words, we are considering a square- netted 
membrane which was brought into thorough discussion in Ref. |[L0|| . Readers may find 
convincing arguments why J(x, y) is a physically sensible statistical measure. In particular, 
the authors think of a polymer chain, whose natural statistical measure is the angles between 
adjacent links. Likewise, for a fluid membrane, they found that the curvatures are the 
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right statistical measure, continuing the polymer into a new space dimension to build up a 
membrane. 

Our theory has the translational invariance of — > + Ah, and thus, the absolute 
value of should be meaningless. Therefore, we introduce the link variable, 

s — adh, (5) 

denoting the "step" at each link; see Fig. [1]. Obviously, we are just performing the well- 



known duality transformation ||24j| , which is very successful in the study of random surfaces. 
Note that now, we arrive at a dual model with step variables {%}; see Fig. |I|. For this dual 
model, the bending-energy cost exists at each plaquette, because it was originally a vertex 
possessing a curvature spanned by adjacent links. Hence, the statistical weight is associated 
at each plaquette; 



with, 



, S2, S3, S4J 



(l + (^) 2 + ( 



^2+^4 ^2 
2a / 



1.5 

(7) 

See also Fig. |I]for the notation of {s Q }. We have set k B T = 1, because k B T can be absorbed 
into the bending rigidity k. It is apparent from the construction of the dual theory that the 
step variables are not fully independent. In the notation of Fig. [I], there exists a restriction 
of si + S2 — S3 — S4 = for each open plaquette. Therefore, we introduce the following 
statistical weight for each of them; 

A(s 1; s 2 , s 3 , s 4 ) = <5 Sl+S2 _ S3 _ S4i o- (8) 

As a consequence, we reach the lattice field theory with the statistical weights p (0) and 
A (|8|), which are arranged in the checker-board pattern. Likewise, the transfer matrix is 
constructed as a strip-like segment shown in Fig. The transfer matrix is diagonalized by 
the density-matrix renormalization group. We will explain it in the next subsection. 

In order to implement the above theory to computer simulation, we must carry out yet 
another simplification. Namely, we discretize the link variable as follows; 

Si = 5 s {i-N s /2- 0.5), (9) 



with i = 1, • • • , N s . The unit of step 5 S is determined selfconsistently through the simulation; 



5 S = Ry/ffl, (10) 

where (• • •) denotes the thermal average. We made several trials for the tuning parameters 
of N s and R; this discretization is the most influential factor concerning the reliability of 
our simulation. We will explore the reliability in the next subsection. 

Our theory resembles the so-called solid-on-solid model, which exhibits Kosterlitz- 
Thouless critical phase. Note that in our model, there is no surface-tension term, and 
there exists the bending elasticity k instead. Therefore, our theory is not right at critical 
phase, but it is rather driven off from it. The renormalization-group flow is actually the 
central concern of the present paper, and it is explored in the next section. 

In the above, we did not pay any attention to the measure factor for the partition sum. As 
is emphasized in Introduction, the measure factor should alter the underlying physics even 
qualitatively. As is apparent from the above formalism, particularly from Eq. (|9]), we accept 
uniform measure for the step variable; namely, we are accepting the normal displacement as 
the statistical measure. Following the idea advocated in Refs. || |10[, we will also consider 
the local mean curvature as for the statistical measure. The replacement of the integration 
variables is absorbed into the redefinition of the statistical weight. Namely, we made the 
replacement, 

dJ{si, s 2 , s 3 , s 4 ) 



p(si, s 2 , s 3 , s 4 ) -> p(si, S 2 , S3, s 4 ) 



\ a=l 



ds a 

The square root is intended to take geometrical mean, because each step variable s a is 
sheared by an adjacent plaquette as well. 

In addition, we consider the local tilt angle as for the measure; namely, we made the 
replacement, 

4 s a 

p(si, s 2 , s 3 , s 4 ) — ► p(si, s 2 , s 3 , s 4 ) n cosatan— . (12) 

a=l a 

As a consequence, we have prepared three types of statistical weights, h (0), J ( p|) and 9 
(|12|), which are examined in the next section. 
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B. Diagonalization of the transfer matrix with density-matrix renormalization 
group 



In the previous section, we have set up the transfer matrix. In principle, one would extract 
various informations, if one could diagonalize the transfer matrix. However, in practice, the 
matrix size exceeds the limit of computer memory size. Such difficulty arises in common 
with such systems called soft matters, which exhibit, by nature, vast numbers of vibration 
modes. Therefore, the Monte-Carlo technique has been employed in order to simulate soft 
matters. However, after the advent of the density-matrix renormalization group |Tl~l , 0) HU > 
the difficulty was removed, and now, soft matters (elastic systems) have come under through 
investigations by means of the diagonalization method; the examples are lattice vibrations 



15, 16], collection of oscillators as a heat bath [17], string meandering motions [IS, 19 



and lattice scalar field theory [p0[| . In essence, the technique allows us to discard "irrelevant 
states," and hence, the number of states are truncated so as to save the computer-memory 
space. 

Below, we will explain the density-matrix renormalization group. Our algorithm is stan- 
dard, and pedagogical guide would be found in the proceeding |H| as well. The transfer 
matrix is represented by a strip shown in Fig. ffl. Our goal is to diagonalize the transfer ma- 
trix with sufficient length. We will show that this goal is achieved by the recursive application 
of the density-matrix renormalization group. The renormalization procedures are presented 
in Fig. |2|, where two sets of renormalizations are shown. Through each renormalization, a 
"block" and the adjacent site are renormalized into a new "block';" and similarly, block' — > 
block". The crucial point is that the number of states for block is kept within m. Therefore, 
we can iterate the procedure until the strip becomes sufficiently long. Such truncation of 
block states within m is managed in the following manner |TT|, |T2||: One first constructs the 



density matrix Pb+« [Q, 0, [13], [14]] with respect to the part of system composed of the block 
and the adjacent site. Then, "relevant states" are extracted from the eigenstates {u a } of 
the density matrix with relatively larger eigenvalues {w a }; the density-matrix eigenvalues 
{w a } are called "weight." Because the weight w a becomes almost negligible for large a, we 
just remain relevant m states with dominant weight, and discard the other remaining states. 
In this manner, the block and the site are renormalized into a new block, whose states are 
represented by the truncated bases {u a } with a — 1, • • • , m. 
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Let us demonstrate the reliability. We will accept the local curvature as the statistical 
measure. In Fig. |3|, we present the distribution of weight {w a } after 40 renormalizations for 
bending rigidity k = 0.5. Simulation parameters are set to be m = 13, N s = 9 and R = 0.9; 
see Eq. (|S]) for details. We see that the weight vanishes very rapidly. In the simulation, we 
retain m = 13 weighted states, and we attain the precision of ~ 10~ 5 . In Fig. we present 
the probability distribution of the step variable. Because we have bounded the range of 
step variable as in Eq. (|9|), we must check whether it covers the actual thermal fluctuation 
deviation. The range of step variables seems to cover the actual fluctuation deviation, hence, 
our treatment of Eq. (f|) turns out to be justified. In all simulations presented hereafter, we 
had monitored such performance check carefully. 



III. REAL-SPACE DECIMATION AND EFFECTIVE BENDING RIGIDITY: AP- 
PLICATION OF THE DENSITY-MATRIX RENORMALIZATION GROUP 

In this section, we study the effective bending rigidity of a fluid membrane. For that 
purpose, we introduce a scheme of real-space decimation p5| . Then, we apply the density- 



matrix renormalization group. All data are calculated after 40 renormalizations; namely, 
the length of the transfer matrix extends to L — 80. 

To avoid confusions, we will comment a few words: In this paper, the word "renormaliza- 
tion" is used in two different contexts. First, we employ the density-matrix renormalization 
group as a simulation technique. Secondly, we manage the real-space renormalization to get 
information of effective bending rigidity. The former terminology is named after the fact 
that we renormalize irrelevant states in order to save computer-memory size. The latter is 
aimed to gain the renormalization-group flow with length scale changed. 



A. Real-space decimation 

Real-space decimation (coarse graining) was first introduced in the study of critical phe- 



nomena 



26]; in particular, it was aimed to interpret the scaling hypothesis. Later on, the 
idea was extended to meet more practical purposes such as quantitative estimation of critical 
exponents and flow equations p7|. We follow Swendsen's version of real-space decimation, 



which has been proven to be very successful p 
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In Fig. ||, we have presented the real-space-decimation procedure. Note that two unit cells 
are renormalized into one enlarged unit cell. In other words, two molecules are renormalized 
into a new molecule, and hence, the parameter M in Eq. (0) is M — 2. Correspond- 
ingly, from microscopic step variables {s a } and {t a }, we construct the coarse-grained step 
variables; 

51 = {si + h)/2 

5 2 = s 3 /2 + s 2 

5 3 = s 3 /2 + t 3 

5 4 = (s 4 + U)/2. (13) 
For this coarse-grained length scale, the curvature is given by, 

( S4-S1 1 S2-S3 \ f-i 1 ( S!+S 4 \2 , ( S 2 +S 4 \2\ _ ( S!+S 4 \2 Si-Si _ ( S 2 +S 4 \2 S 2 -S 3 
j_ l(L5ap 3 (1.5a) 2 H 1 I 3a ' t I 3a ' ) \ 3a ' (1.5a) 2 ^ 3a ' (1.5a) 2 

After that coarse-graining, one must redefine the unit of length accordingly; namely, 

a — > a/ \pl 

hij -> ^7^2. (15) 

The coarse-grained membrane may be described by the Hamiltonian H = J dAn'J 2 /2 
with renormalized bending rigidity k' . As is well known, the transformation coefficient 
dn'/dn contains much information on the infrared behavior of effective coupling constant 
p5| . Anticipated behaviors of Kj— fx £1X6 drawn schematically in Fig. |6|. From the figure, 
we see that for dn'/dn > 1, the membrane is stiffened in the infrared limit, whereas for 
dn'/dn > 1, the membrane is softened. The transformation coefficient Ok'/Ok is given by 



the chain relation 25 



d(J 2 dA) _ dK' d(j 2 dA) 

8k Ok 8k' 

Here, dA denotes the area of the membrane segment shown in Fig. H Remaining task is 
to perform the above numerical derivatives. Numerical differentiations are possible, because 
our data are free from statistical errors. That is a great advantage of our algorithm over 
others such as Monte Carlo. We adopted the "Richardson's deferred approach to the limit" 
algorithm in the text book 28|| . In this algorithm, one takes an extrapolation after calculat- 



ing various finite-difference differentiations. We monitored the relative error, and checked 
that the error is kept within 10~ 3 . 



B. Membrane stiffening in the case of local-curvature measure 



Following the idea of Refs. || |10[, we will accept the local curvature J(x,y) as for 



the statistical measure. This is achieved by adopting the statistical weight (|ll|) in the 
construction of transfer matrix. In Fig. [?|, we plotted the transformation coefficient Ok' /dn 
(p"6|) for the bare coupling 0.17 < k < 4. In obtaining the data, we have made a number 
of trials for the simulation parameters of m (number of block states), N s and R (range of 
step variables); see Eq. (^|) for details. Such technical informations are also presented in the 
figure caption. Thereby, we confirm that good convergence is achieved with respect to the 
tuning parameters. 

From Fig. [7], we see that the inequality 8k' /dn > 1 holds. Referring to the anticipated 
behavior depicted in Fig. |6], we found that the membrane is stiffened effectively for macro- 
scopic length scales. As a matter of fact, for other integration measures, the inequality 
dn' j <9k > 1 is not satisfied as would be shown in the next subsection. Our simulation result 



is the first ah initio support of the analytical argument by Pinnow and Helfrich H |10 1 . 

In the figure, for small rigidity k < 1, where thermal undulations should be enhanced, we 
observe notable enhancement of the transformation coefficient cW /dn. Hence, it is indicated 
that this stiffening is driven, quite contrary to our naive expectation, by thermal undulations. 
This rather counterintuitive result suggests that the "hat excitation" ||, which is a sort of 
solitonic excitation, would be created over the membrane surface: The hat excitation is a 
solitonic object accompanying localized dimple-like deformation. In Ref. 0, the author 
claimed that the hat excitations give rise to, unlike conventional sinusoidal undulations, 
the membrane stiffening. It would be astonishing that for such small rigidity k < 0.5, 
where the membrane should be crumpled considerably, the concept of hat excitation is still 
applicable. This fact may tell that the hat excitation is indeed solitonic in the sense that the 
hat excitations are stable under collisions, and the thermal undulations can be decomposed 
into individually propagating hat excitations. 

It is to be noted that such pronounced enhancement of dn'/dn for small k is not cap- 
tured by the analytical one-loop renormalization-group treatment (Q), because it just yields 
dn'/dn = 1 for over all v /t. (It might be convincing, because analytical perturbative treat- 
ment should be justified for sufficiently rigid membranes.) This is obviously the advantage 
of our first-principle simulation. For exceedingly small rigidity k < 0.17, however, the mem- 
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brane becomes too much crumpled, and we cannot continue reliable simulation. (From a 
technical viewpoint, the density- matrix weight exhibits almost uniform distribution, and we 
cannot set any reasonable threshold for the truncation of states.) 

On the other hand, for large-rigidity side n > 1, we see that the transformation coefficient 
approach the neutral value Dk' /dn ~ 1. Note that the analytical one- loop renormalizat ion- 
group analysis (0) yields dn'/dn = 1. Hence, we see that for sufficiently large k, the 
analytical result (H) is recovered asymptotically. 



C. Membrane softening under the statistical measures of h(x,y) and 9{x,y) 

For other integration measures such as normal displacement h(x, y) (0) and tilt angle 
6(x,y) (0), we found behaviors quite contrastive to that shown in the previous subsec- 
tion. In Fig. [8], we plotted transformation coefficient dn'/dn for the normal-displacement 
integration measure (0). (Note that this integration measure has been widely used so far 



in the analytical calculations except Refs. H, ITOf.) From the figure, we see definitely that 



dn 1 /dn < 1 holds. Hence, in this case, the membrane is softened effectively for macroscopic 
length scales; see the schematic behavior shown in Fig. |6|. In addition, it is to be noted 
that for small rigidity k < 1, the transformation coefficient dn'/dn is suppressed, and it 
approaches the neutral value dn'/dn = 1 as the membrane rigidity k increases. This fact 
tells that the softening is driven by the thermal undulations. Note that the effect of ther- 
mal undulations appears in a way quite opposite to the aforementioned J(x, y) statistical 
measure. 



In the case of tilt-angle statistical measure (|12D, the transformation coefficient exhibits 
similar behavior dn'/dn < 1; see Fig. |]. However, Ok' /dn is much closer to the neutral 
value dn' /dn ph 1, suggesting that the extent of softening is smaller than that of the normal- 
displacement integration measure. In other words, the membrane shape would stay almost 
scale-invariant under the choice of tilt angle as for the statistical measure. 

In the calculations presented in this subsection (particularly, for the h(x,y) measure), 
the simulations suffered from instabilities coming from unbounded undulations due to the 
membrane softening. More specifically, during the simulation, membrane becomes crumpled 
spontaneously, and the membrane shape is trapped to a certain metastable configuration. 
Such pathology may arise, because there are exceedingly numerous thermally activated 
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configurations of almost equal statistical weight. Diagonalization of transfer matrix thus 
fails in searching true globally-stable thermal equilibrium. Moreover, the postulation of Eq. 
(|) would not be fully justified, where we had assumed that the range of step variable is 
bounded. As the simulation parameters of m and N s are improved, those instabilities are 
avoided to some extent. In a sense, those instabilities reflect the fact that the membrane is 
softened indeed, and the shape fluctuations are enhanced for long distances. From the above 
experience, we are led to the speculation that the Monge gauge would not be very justified 
for the description of fluid membranes in the case of softening, at least, beyond perturbative 
level. 

IV. SUMMARY AND DISCUSSIONS 

We have investigated the effective bending rigidity (|2]) of a fluid membrane ([[]) for macro- 
scopic length scales. The effective rigidity has been arousing renewed interest, since Pinnow 
and Helfrich pointed out that the membranes would be stiffened by thermal undulations. 
The key ingredient of their argument is that the local curvature should be the right statisti- 
cal measure for the partition sum. Motivated by this remarkable scenario, we had performed 
first-principle simulation with the transfer-matrix method and the density-matrix renormal- 
ization group; see Figs. [I] and |2|. Our simulation scheme does not rely on any perturbative 
treatments, and it covers various statistical measures such as local curvature (|TT|) , normal 
displacement (0) and tilt angle (|l^) . Because analytical variable replacements among those 
variables require rather tedious task even on perturbative level, it would be meaningful to 
survey various measure factors systematically by ab initio simulation. Performance of the 
simulation scheme is demonstrated in Figs. |3] and |]. 

We introduced a scheme of real-space decimation as is shown in Fig. [5], and correspond- 
ingly, we defined coarse-grained curvature J (|H]). Those preparations enable us to calculate 
the transformation coefficient dn'/dn ([16]), from which we read off the direction of renor- 
malization group flow; its anticipated behaviors are drawn schematically in Fig. As a 
consequence, in respect of statistical measures, we observed clear distinction between the 
mean curvature dn' /dn > 1 and the others dn' /dn < 1. Namely, when we accept mean cur- 
vature as a statistical measure, the effective rigidity is, in fact, renormalized toward strong 
coupling; see Fig. [?]. That is, membrane stiffening takes place. On the contrary, both normal 
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displacement and tilt angle appear to lead to membrane softening; see Figs. || and |9|. Our 
simulation results are the first ab initio support of the aforementioned analytical argument 

The membrane stiffening, apart from its mere curiousness, is quite favorable for the con- 
sistency of our numerical simulation: In our simulation, in the first place, we had restricted 
the range of step variables as in Eq. (]T5|) . Apparently, this postulation is in favor of mem- 
brane stiffening. In addition, we had discarded irrelevant states through the density-matrix 
renormalization group in order to keep the number of states tractable in computers. Again, 
this truncation of states is validated consequently after the onset of membrane stiffening. 
Furthermore, we believe that the very starting point of our theory, namely, the Monge gauge, 
is in favor of membrane stiffening, because it is assumed that there exists a reference plane 
from which all undulations are excited. Membrane stiffening may serve solid grounds for 
the Monge gauge. 

For those reasons, we could perform well-controlled simulation in the case of membrane 
stiffening. On the contrary, as for the case of membrane softening, our simulation faces a 
number of conflicts coming from unbounded undulations as the strip length of transfer matrix 
is enlarged; see Section |111 C| for details. Because such pronounced undulations deprive the 
grounds of the reference plane, it would be rather suspicious to adopt the Monge gauge 
at least beyond perturbative level. By the way, we are fairly confident of the efficiency 
of our simulation in the case of membrane stiffening, and we believe that the membrane 
stiffening may serve a new promising research field for the application of the density-matrix 
renormalization group. 

The following problems remain open: First, In the article f8|, the author raised an in- 
triguing picture of deformation mode called "hat excitation." This is a sort of soliton which 
is in contrast to the naive sinusoidal undulation. With the hat-excitation picture, the author 
gave valuable physical insights. First-principle examination of the reality of this hat excita- 
tion would be desirable. Secondly, irrespective of the membrane stiffening, the correlation 
length is known to be finite, and the orientational correlation is lost for long distances j9]. 
In Ref. 0, it is speculated that those two modes, namely, orientation and curvature, would 
be decoupled. A deeper understanding of the above points would establish the justification 
that the local curvature is the right statistical measure. 
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(c) 

FIG. 1: (a) On the square lattice, we consider scalar field hij denoting normal displacement of a 
membrane with respect to a reference plane. Step variable (gradient field) s = adh is defined at 
each link, (b) The local statistical weights p (R) and A (||) are represented by shaded and open 
squares, respectively. The statistical weight p has several variants so as to take account of other 
integration measures such as curvature (|Tl| ) and tilt angle (|l^). (c) From those local statistical 
weights, we construct a strip whose row-to-row statistical weight yields the transfer- matrix element. 
This transfer matrix is subjected to the DMRG diagonalization as is shown in Fig. 0. 



16 



block 



DMRG 



block' 



y 
block' 



block" 



block 



block' 



DMRG 



» block"->block 



FIG. 2: Schematic drawing of the density-matrix normalization group (DMRG) procedure [11 



12, 13, 14]. From the drawing, we see that through DMRG, a "block" and the adjacent site is 
renormalized into a new "block"'; similarly, bUdjck' — > block". At this time, the number of block 
states is retained within m; see text. In this manner, we can diagonalize large-scale transfer matrix 
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FIG. 3: Distribution of the eigenvalues (weights) {w a } (a: integer index) of the density matrix for 
bending rigidity k = 0.5. Simulation parameters are m = 13, N s = 9 and R = 0.9; see Eq. @. 
We see that w a vanishes very rapidly for large a. In the simulation, we remained relevant states 
up to m = 13, and we achieved the precision of ~ 10~ 5 . In this manner, the number of states of 
"block" in Fig. ^ is truncated (renormalized). 
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FIG. 4: Probability distribution of the step variable P(si); see Eq. (Q), where we had set threshold 
for the range of step variable like {sj} (i = 1, • • • , N s ). The simulation parameters are the same as 
those of Fig. § The range of step variables seems to cover the actual fluctuation deviation. 
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FIG. 6: Schematic flow diagram of the effective bending rigidity k — * k'. Depending on the 
transformation coefficient 8k' /8k > 1 (< 1), the membrane is stiffened (softened) in the infrared 
limit. 
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FIG. 7: Transformation coefficient 8k' /dn is plotted for bare rigidity k. We have accepted the 
local curvature as for the statistical measure. The simulation parameters for each symbol are (+) 
m = 12, N s = 7 and R = 0.9, (x) m = 12, N s = 7 and R = 1.1 (*) m = 11, N s = 8 and ft = 0.9 
and (□) m = 15, iV s = 6 and i? = 1.1; see Eq. (^) for detail. Referring to the anticipated behavior 
shown in Fig. |(| we see that the membrane is stiffened effectively in the infrared limit. 
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FIG. 8: Transformation coefficient dn'/dn is plotted for bare rigidity k. We have accepted the 
normal displacement as for the statistical measure. The simulation parameters for each symbol are 
(+) m = 12, N s = 7 and R = 0.8, (x) m = 10, N s = 9 and R = 0.6 and (*) m = 11, iV s = 8 and 
i? = 0.7; see Eq. (|9|) for detail. Referring to the anticipated behavior shown in Fig. ||, we see that 
the membrane is softened effectively in the infrared limit. 
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FIG. 9: Transformation coefficient dn'/dn is plotted for bare rigidity k. We have accepted the 
local tilt angle as for the statistical measure. The simulation parameters for each symbol are (+) 
m = 10, N s = 9 and R = 0.7, (x) m = 12, N s = 7 and R = 0.7 and (*) m = 11, N s = 8 and 
R = 0.7; see Eq. (|9|) for detail. Referring to the anticipated behavior shown in Fig. ||, we see that 
the membrane is softened effectively in the infrared limit. 
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